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Abstract[] On the basis of Watt’ s mathematical model of the index of population trend] /[] an approximate variance 
estimator of J was developed with delta-method in this paper. Taking the natural population life tables of rice leaf 
roller Cnaphalocrocis medinalis{] Lepidopteral] Pyralidae[] published by Wu et al .[ 1986[] as an example[] it was 
applied in evaluating the control effectiveness of biological agen{] Trichogramma japonicum Ashmead[] and chemical 
insecticidd] powder{] 1.5% mevinphos + 3% alpha-hexachloro cyclohexane]. By Z-test criterion[] biologically and 
statistically sound conclusions were drawn that at the generation level[] the suppression effectiveness of trichogramma 
wasp against rice leaf roller was better than the blank control with P =0.0111 and smaller 7 valud] 0.0390 versus 
0.1768L[] and it was also statistically better than the insecticides with P = 0.0036 and smaller J valued] 0.0390 


versus 0.3035 

















. The insecticides was not statistically worse than the blank control with greater 7 value] 0.3035 


versus 0.1768L[] and P = 0.2236. Although these conclusions are close to that drawn by Wu et al.) 1986[] only 
based on the index of population trend point estimates[] they are supported by an accuracy measurement P-value. 
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medinalis [|] delta- method 


Since Morris and Miller{] 1954[] introduced life 
table methodology into insect natural population 
dynamics analysis[] this method has been widely applied 
in evaluating the effectiveness of physical or biological 
factors on insect population fluctuation[] and also in 
determining the management effectiveness of chemical or 
biological agents in agricultural pest management 
practice. Watt] 1961[] 1963[] developed mathematical 
models computing the index of population trend] 7/00 
which provided an integrated index containing all 
regulation effect information of environmental factors on 
the whole generation of the target population. Harcout 

O 1963[] 1969[] discussed life table construction tactics] 
its application scenario[] and future applications in 
details. Pang et al [] 1981[ discussed field experiment 
design[] data collection[] and survival rate estimation 
tactics for constructing a type of natural population life 
table of rice leaf roller] Cnaphalocrocis medinalis [] 
and its application in determining the effectiveness of 
environmental factors. By the rice leaf roller natural 
population life tables[] Pang et al .[] 1984[] used index 
of population trend as the analytical tool evaluating the 
management effectiveness of chemical and biological 
agents in its IPM system at generation level] and 





concluded that at generation level biological agents were 
not only environmental friendly but also more effective 
and efficient than chemical insecticides. Bellows et al . 
O 1992[] and Carey[] 2001 [] reviewed the recent 
development in life table methodology focusing on 
experimental design[] data collection[] and mortality 
estimation[] and its application in analyzing insect 
population dynamics[] estimating the effectiveness of 
environmental factors[] such as natural enemies and host 
plant traits[] and determining control effectiveness of 
biological and physical factors in pest management. 
Unfortunatelyl] up to date there is still no publication 
available focusing on the variance estimation of index of 
population trend] /[]. By Watt’ s mathematical model 
and delta-method[] the approximate variance estimation 
of the index of population trend and its application are 
the topics concerned in the following sections. 


1 AN APPROXIMATE VARIANCE ESTIMATOR 
FOR INDEX OF POPULATION TREND] J[] 


1.1 Mathematical model of index of population 
trend] Z0 
By the mathematical model developed by Watt 
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O 19610 1963[[] Pang et al [] 1984[] applied the 
following equation computing the index of population 
trend] /{] for natural population life table of rice leaf 
roller[] 

T= S,S,AS,S,5/S,N. 
Where S, is the survival rate at egg stage[] S,..-[]S,, is 
first ...1] last stage 
respectively[] S, is the survival rate at pupa stage[] S; is 


the survival rate at larva 
the female moth rate[] S, is the proportion of female 
moth achieving standard egg number calculated with the 
actual mean egg number divided by theoretical standard 
egg numberl] N, is the standard egg number per female 
moth. 
1.2 Delta-method 

Delta-method based on Taylor expansion is often 
applied to estimate the approximate variance and the 
corresponding confidence intervals for a random variable 
for which there is not explicitly mathematical formula 
available to estimate its variance [| Sorensen and 
Gianola[] 2002[]. Let x[] y be two random variables 
with the following function relation[] 

y = As 

And assume that this function is differentiable at least 
twice. Let xọ be the expected value of x. If the value 
of F] x[] is known at the point of xoL] denoted by 
[A x LL] the value of y at aL] which is pretty close to its 
expectant[] can be approximately formulated by Taylor 
expansion as follows. 


fa xO =f xo + f xol x — xol 
+30 xoLL] x- xo 


Remove the constant term [f] xo[ to the left of the equal 
sign of the previous equation[] there is 


A x0- xo = xoL] x — xol 
+50 xof x — x 


If the variance of random variable x is known[] the 
variance of random variable y can be approximately 
estimated by the following mathematical formula through 
only retaining the second order term[] 
vad yO = vl A <O-A DO A xo vdd «OO 20 
Or if the variance of random variable y is known[] the 
variance of random variable x can be approximately 
estimated by 
vad xO <0 A xf vd yO «.... O 30 
1.3 An approximate variance estimator for index 
of population trend] Z[] 
Transforming equation [| 1[] of the index of 
population trend] /{] with logarithm[] there is 
lẹ 10 = ldd S.O + lob) S,\0+ A + log 5,0 
+ log S,0 + ldg S/O + log 5,0 
Filda Wali eam seats bao 040 


By delta-method[] the variance of log-transformed 
survival rate is easily calculated[] 

vdi log SiO ~ vd SDS; 
Since the survival rate of released cohort of rice leaf 
roller is related to binomial distribution[] the variance of 
the survival rate] S;[] is computed with the following 
formula[] 

asg 21-50 

nj 

Here n; is the number of individuals released at ith life 
stage. Actually[] the true survival rate is unknown. The 
maximum likelihood estimate of survival rate[] $; [0 is 
substituted into previous mathematical formula[] and the 
maximum likelihood variance estimate of the survival 
rate is gotten. With similar procedure as survival rate[] 
the variance of female moth rate is estimated. 

For the reproductive component at adult stage[] 
denote the mean egg number per female moth and the 
corresponding variance as x[ and vat] x00 which are 
obtained from cohort experiment. The proportion of 


female moth achieving standard egg number is S, = = 


N.’ 
vad x0 
~ The 


transformed female moth 





The corresponding variance is vaf] S, O= 


variance of logarithm 


proportion achieving standard egg number is 
Ney vad x0 

vd Idd S.Q = (=) vail S.O= = 

x 


X 





By the field experiment design and data collection 
tactics with 
independence among sequential life stages is assumed 
for simplicity. Thus[ 


aidie o Sa ee So 


Jei Tel 
Where m is the total life stages[] n; is components in 


release and recapture method[] the 


jth life stage. The standard egg number for a female 
moth] N, [O a constant[] vanishes in the variance 
estimation model. By the variance estimates for ldg ZOO 
the variance of the index of population trend] Z[] is 
vad JQ = exp log IOF vdd ldg 700 
= Pvd log 10 
Since the true value of index of population trend is also 
unknown[] after substituting its maximum likelihood 
estimatd] Ê[ into this equation[] its maximum likelihood 
variance estimator is[] 
va] i] = exp log] fvdd lol 70 
= Pv log 10 
1.4 Confidence intervals for index of population 
trend] Z[] and Z-test criterion 
By central limit theorem[] when the sample size is 
large enough[] the following statistics approximately 
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tends to standard normal distribution 
fod n>> 
> M oo 
Jv vay 10 
Where n is the sample size[] J the true value of index of 
population trend[] and N] O[J1{] the standard normal 
distribution. By this law[] the 1 — @ confidence intervals 





of the index of population trend] Z[] is estimated as 


Oi- Z, xv vdi (001 + Z, xv vad 700 ... 
Here Z 7 is the critical value from standard normal 


Puun a 
distribution at 1 - >. 


2 
Depending on the null hypothesis that there is no 
difference of index of population trend between two 
treatments[] J; - I = OCU] and the independence 
between treatments guaranteed by the experimental 
design{] the following Z-test criterion is applied in 


treatment comparison[] i 
L—L 
z-— a iU =F eaaa] O 100 
vd] 40 + vd 20 
By this Z-value[ for two tail test[ the corresponding P- 
value is easily obtained. 


2 APPLICATION OF THE APPROXIMATE 
VARIANCE ESTIMATOR 


2.1 Life stage specific life table of rice leaf roller 
and the hypothesis test 

Based on the life tables of rice leaf roller from 
three field treatments[] blank control plot biological 








control plot released with Trichogramma japonicum 
Ashmead[] and insecticide control plot dusted with 
1.5% mevinphos + 3% alpha-hexachloro cyclohexane 
powder[] published by Wu et al.[] 1986[]] their 
stage table[] the 


maximum likelihood estimate of the index of population 


corresponding life specific life 
trend] Z[] for each treatment[] and its approximate 
variance estimate are listed in Table 1. In insect 
population management experiment|] the null hypothesis 
O H,O assumes that the control agent does not change 
survival and reproduction patterns of the target host 
population or the control agent does not change the 
index of population trend] Z[[] which is tested with life 
table methodology. The alternative hypothesi{] H, U is 
that the control agent significantly changes the survival 
and reproduction patterns of the target host population or 
it significantly changes the index of population trend 
with type I error probability a = 0.20 0.10 0.050 or 
0.01 respectively. By Z-test criterion[] the index of 
population trend] /[] in each treatment and their 80% [] 
90%] 95 %[] 99% confidence intervals are plotted in 
Fig. 1. By the confidence intervals of Z[] when the type 
I error a@ = 0.05[] the biological control plot was 


statistically better than the blank control and the 
insecticide treated plot. When the type I error a = 
0.01[] there was no difference among all the three 
treatments. By the comparison between treatments with 
Z-test criterion[] the suppression effect of releasing 
wasp was statistically better than blank control with 
smaller [-valud] 0.0390 versus 0.1768[] and P-value 
O 0.01110 < 0.05[] and also statistically better the 
chemical insecticide with smaller J-value[] 0.0390 
versus 0.3035[] and P-valud] 0.00360 < 0.01. The 
chemical insecticide was not statistically worse than 
blank control with greater J-valud] 0.3035[1] and P- 
valud] 0.2236[] > 0.05. Between these to procedures[] 
the confidence interval comparison is much more 
conservative. Although such results are close to that 
drawn only by the point estimates of index of population 
trend[] these conclusions are supported by an accuracy 
measurement P-value. 


Table 1 Natural population life table of rice leaf roller 
Life Stage 


or Parameter 


Survival rate 





Blank control Trichogramma wasp Insecticides 





Egg 0.4130 0.1850 0.4160 

lst instar 0.4400 0.3400 0.2000 

2nd instar 0.5000 0.3100 0.6200 

Larva 3rd instar 0.3000 0.3300 0.5400 

4th instar 0.7000 0.5600 0.7600 

Sth instar 0.5500 0.6000 0.6500 

Pupa 0.4100 0.3800 0.5400 

Sp 0.5000 0.5000 0.5000 

Adult S, 0.4100 0.4100 0.4100 
Ne 200 200 200 

i 0.1768 0.0390 0.3035 

Variance of 7 0.0028 0.0002 0.0082 


3 CONCLUSIONS 


In data analysis the point estimates for 
population parameters are useful to show population 
properties|] such as mean[] median[] quartiles[] ete. 
When such estimates are applied to compare properties 
of different populations[] such as insect on different host 
plants or on different resistance level cultivars[] in 
different geographical areas[] and with different genetics 
traits[] it is impossible to draw statistically dependable 
conclusions only by point estimates[] although such a 
comparison does provide certain information. Thus[] 
how to estimate the variance corresponding to the point 
estimate is a pragmatic problem related to management 
agent selection in pest IPM system. Unfortunately[] for 
some statistics[] like the index of population trend] 700 
it is very difficult to find an explicit mathematical model 
to calculate its variance[] up to date there is still no 
such estimator available. Practically[] it is desirable 
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Fig. 1 


using approximate estimator achieving variance estimates 
for population point estimates. By the analysis and 
discussion in foregoing sections[] the following 
conclusions are drawnl] 

O 10 After ignoring the correlations between the 
sequential life stages of the target insect population by 
the field experiment design{] sampling schedule[] and 
data collection tactics[] an approximate variance 
estimator for the index of population trend is obtained by 
the log-transformed Watt’ s mathematical model and 
delta-method. 

O 20 By the central limit lawl] when the sample size 
is large enough[] the estimates of the index of population 
trend] /{] tends to confirm with M 100° [0 where o° is 
variance estimated by the 
estimator. By this law[[] 1 — a[] confidence intervals 
can be estimated for each point estimate of the index of 


population trend[] and the Z-test criterion can be 


approximate variance 


applied into the evaluation of the pest management 
agents. 

O 30 By Z-test criterion[] the biological control plot 
was statistically superior to the blank control] P = 
0.0111[] and the chemical treated plo{] P = 0.003600 
but the chemical treated plot did not differ significantly 
0.2236. Although 


from control with P = such 


oor alpha=0.1 
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Index of population trend and its 1 — a confidence intervals in control[] wasp] and insecticide treated plots 


conclusions are close to that drawn by Wu et al. 
O 1986[] only based on the point estimates of index of 
population trend[] they are supported by an accuracy 
measurement P-value and statistically dependable . 

This approximate variance estimator for index of 
population trend{] Z[] is based on the following 
assumptiong ||] a[] the distribution of the index of 
population trend] /[] tends to normal when sample size 
is large enough[] and[] b[] the correlation between 
survival rates of the sequential life stages is negligible. 
Actually] it is very difficult to answer whether the 
normality of the estimates is achieved and the 
correlations between the survival rates of sequential life 
stages are negligible. By biological relations within the 
interaction system composed of rice plant[] rice leaf 
roller] its natural enemy[] and its food competitors[] 
more or less there exists direct or indirect correlation 
between survival rates of the sequential life stages. 
Thus 
techniques for the index of population trend] Z[ based 
on life table data is still a question worthy of future 
investigation . 


searching for rational variance estimating 
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